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In this article we report on a new proposal to treat the infrared problems of thermal QCD by 
taking into account explicitly the confining influence of the Gribov horizon. In order to make 
clear the possible value of such an approach, we briefly review the most important arguments why 
a straightforward perturbative description of finite-temperature QCD is unlikely to be successful. 
From the infrared problems of thermal perturbation theory one can conclude that confinement 
effects and bound states probably play an important role also in the high-temperature phase. 
To set the stage we recount the supposed role of the Gribov horizon for confinement, before we 
turn to the application to finite-temperature theory. In the current approach it has been found that 
the contributions to the free energy from the explicit inclusion of the horizon begin to set in at 
order g 6 - precisely where the infrared problems of thermal QCD lead to a breakdown of ordinary 
perturbation theory. 

From the study of observables (free energy, anomaly, bulk viscosity) we also note that for ther- 
modynamic observables the leading order term obtained by such an expansion in the coupling 
strongly deviates from the more complete numerical solution. This can be regarded as yet another 
sign for general problems of series expansions in thermal QCD. 
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1. Infrared Problems in Thermal QCD 

The properties of QCD at high temperatures (^(100 MeV) and higher) are obviously impor- 
tant for an understanding of the early universe. While this domain is to some extent accessible 
experimentally (at RHIC and LHC), it remains a challenge for theory. 

Initially one had expected to find a phase transition where at some temperature hadrons melt 
into a "quark-gluon plasma" [Q]. In such a plasma, quarks and gluons were supposed to be almost 
free particles which could be reliably described by perturbation theory. Such a purely perturba- 
tive description was motivated by the fact that high temperatures imply high average momentum 
transfer /J. and thus - due to asymptotic freedom - a small coupling g(ju). 

However, this picture turned out to be too naive. In principle, this should have been clear at 
least since 1980, when it was shown |^, |3]] that at order g 6 a natural barrier arises for any per- 
turbative description (the Linde problem). Even earlier than that, the simple fact that the infinite- 
temperature limit of four-dimensional Yang-Mills theory is a three-dimensional confining Yang- 
Mills theory could and should have been regarded as a sign that any straightforward perturbative 
approach to high-temperature QCD was problematic. 

Still it took more than 20 years until it began to be accepted that the high-temperature phase of 
QCD has little to do with a conventional plasma. The results of the RHIC experiments [Q] clearly 
showed that also above the transition, bound state phenomena can not be neglected. 

Since perturbation theory is limited to a finite order (and, in addition, for experimentally acces- 
sible temperatures the apparant convergence of the perturbation series is bad []|, ^]), other methods 
have been invoked to treat the properties of finite temperature QCD: 

One popular proposal ^ is to explicitly separate the three relevant scales (hard - 2nT, chro- 
moelectric - gT and chromomagnetic - g 2 T) by introducing effective theories at each scale. Two 
of these theories can be treated perturbatively, while one (which governs the magnetostatic sector) 
is genuinely nonperturbative. 

Therefore nonperturbative methods are - either in the direct or in the effective theory approach 
- urgently needed. One possibility is given by functional methods, based for example on Dyson- 
Schwinger equations; for an application to finite temperature see e.g. [Q, |9]]. These methods provide 
valuable insights, but unfortunately up to now the pressure (and quantities derived from it) are 
difficult to access. Still functional methods provide further evidence for the picture of bound states 
playing an important role even at very high temperatures and parts of the gluon spectrum being 
confined at any temperature. 

The nonperturbative methods most actively pursued for gauge theories are lattice-based, and 
indeed lattice data is available for finite temperatures up to T « 5Tc- In this region for pure SU(3) 
gauge theory the problem of determining the equation of state is regarded as solved Also 



the AdS/CFT duality [11] and the AdS/QCD duality [12] are employed to study finite temperature 
gauge theories. 

In this article we report on a project which follows an alternative path: Since the footprints 
of confinement are visible also in the "deconfmed" phase, the confining influence of the Gribov 



horizon is explicitly taken into account. The method has been proposed in [ ]13| , |14j] and further 
developed in Jl5| , [l6| ]. To explain it, we will briefly summarize the importance of the Gribov 
horizon in section ||, while in section ^] we will discuss the application to finite-temperature theory. 
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Figure 1: (a) In non-Abelian gauge theories there are typically configurations A] , k = 0, 1,2,... on the 
same gauge orbit (9, which fulfill the same gauge condition F[A] = 0. (b) The gluon Dyson-Schwinger 
equation (only the infrared-dominant part is displayed on the rhs). 



2. The Gribov Horizon and Confinement 



The initial goal of gauge-fixing was to single out precisely one configuration on each gauge 
orbit in order to embed the physical configuration space into the space of all field configurations. 



As it has been shown by Gribov [ Jl7[ | and further discussed by Singer [[18], in non-Abelian gauge 
theories the usual gauge-fixing procedure does not yield uniqueness. So for any local gauge-fixing 
condition F[A] = one typically finds, as sketched in figure [j]a, several gauge-equivalent configu- 
rations A(°), AW = ^A(°) := g^A^lgi + g^dgu A( 2 ) = «A(°), . . . which all fulfill this condition. 

This does not affect perturbation theory (since only the vicinity of A = plays a role there), 
but in order to account for nonperturbative effects, one has to consider the effect of such Gribov 
copies. In particular the infrared properties are dramatically changed. 

In order to eliminate Gribov copies one typically restricts the domain of integration to the 
Gribov region, i.e. the region where the Faddeev-Popov operator ^#(A) is positive semidefinite. 1 
The boundary of the Gribov region (where ^#(A) has at least one zero eigenvalue) is called the 
Gribov horizon. The presence of the Gribov horizon (and the restriction of the functional integral 



to the Gribov region) leads to modifications of the Faddeev-Popov procedure [20]. The influence 



of the Gribov horizon can be incorporated in a local action [ 2 1 



While the Gribov problem seems to be a nuisance at first glance, in fact a whole confinement 



scenario is based on it - proposed in [ ]17[ ] and further elaborated by one of the authors [ |21| , |22| ] : 
Since the Gribov region is (on the lattice) a high- or (in the continuum) an infinite-dimensional 
region, due to geometrical reasons most of the volume is concentrated close to the boundary ("the 
entropy argument"). So most configurations which contribute to the path integral have an almost 
vanishing Faddeev-Popov operator and correspondingly an enhanced ghost propagator. 

Such an enhanced ghost can be made responsible for the suppression of the gluon propagator 
in the infrared. This is most easily understood by analyzing the gluon Dyson-Schwinger equation 



(see figure |l|b) [|23|, |24|, \t^]. The ghost loop on the right-hand side is infrared enhanced and thus 



the gluon propagator is accordingly suppressed - the gluon cannot propagate over long distances. 



'However, as discussed in Jl9|], in principle one has to restrict the integration even further - to the Fundamental 
Modular Region (FMR), which is unfortunately characterized by a global condition and thus difficult to handle. 
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(a) 




(b) 




Figure 2: (a) the rescaled Gribov mass m* — j, (b) the rescaled free energy - ^l SB , 



where the Stefan- 



Boltzmann term wsb has been subtracted, [solid - numerical solution, dotted - leading order expansion] 



3. Application to Finite Temperature Theory 



We now turn to the "marriage" of finite temperature theory and Gribov's confinement sce- 
nario. The basic physical idea is that the infrared divergences of finite-temperature perturbation 
theory (which are responsible for the Linde problem) do not arise when the domain of functional 
integration is cut off at the Gribov horizon. The cut-off is done in Coulomb gauge which is well 
adapted to finite-temperature calculations. 

Technically the cut-off at the Gribov horizon is implemented by adding a "horizon function" to 
the action 26]. The initially non-local term then gets replaced by a local, renormalizable term 
in the action by means of an integration over a multiplet of auxiliary ghost fields. The new term 
in the action depends on a mass parameter m; the functional cut-off at the Gribov horizon imposes 
the condition that the free energy W or quantum effective action T be stationary with respect to that 
mass, sJE = -I 1 =0. 

' dm dm 

This "horizon condition" has the form of a non-perturbative gap equation that determines the 
Gribov mass m = m(T \Aqcd) and thereby provides a new vacuum, around which a perturbative 
expansion is again possible. Lowest-order expansion gives the Gribov-type dispersion relation 



Elk 




(3.1) 



From knowledge of the Gribov mass one can deduce contributions to the free energy w, the en- 
ergy e, the (rescaled) anomaly A r = and (making use of a relation obtained in [Q) the bulk 



viscosity £. Results for m and w are displayed in figure ^[ 

An expansion of the results in powers of g is possible at least in leading order, but this leading 
order term strongly deviates from the full numerical result. This is probably not a fault of the 
method, but should be regarded as yet another sign that series expansions in thermal QCD have to 
be treated with great care. 

Nevertheless as a matter of principle, it is a significant success that for thermodynamic ob- 
servables this procedure gives finite results precisely at the order, g 6 at which ordinary perturbation 
theory diverges. It also confirms the picture of thermal QCD composed of two sectors with one 
being accessible perturbatively, while the other one is genuinely nonperturbative. 
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